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Abstract. For the two-dimensional electron gas, the exact high-density limit of the correlation 
energy is evaluated here numerically for all values of the spin polarization. The result is spin- 
resolved into It? Tl) an d || contributions and parametrized analytically. Interaction-strength 
interpolation yields a simple model (LSD) for the correlation energy at finite densities. 

In recent years, two-dimensional (2D) electron systems have become the subject of exten- 
sive research |L. The 2D version of density functional theory (DFT) has proven particu- 
larly successful in studying quantum dots (2JIHJII!- The local spin-density approximation 
(LSD) of DFT requires the correlation energy of the spin-polarized uniform electron gas. 
This quantity in 2D is known accurately for a wide range of densities and spin polar- 
izations from fixed-node diffusion Monte Carlo simulations [Sj. Its high-density limit is 
known exactly in terms of six- dimensional momentum-space integrals Resolved into 
contributions due to ||, ||, and || excitation electron pairs, these integrals are evaluated 
here numerically. The analytical parametrization of the results, Eqs. f!16|) and ()17|) below, 
is a crucial ingredient for the construction of the spin-resolved correlation energy at finite 
densities, performed recently for the 3D electron gas [Zj. It is also required for studying 
the magnetic response of the spin-polarized 2D electron gas [HUH]- Generally, it provides 
a fundamental test for numerical parametrizations of the correlation energy [Sj. 

In the 2D uniform electron gas, the electrons are moving on a plane at uniform density 
p = [7r(r s aB) 2 ] _1 , where 05 = 0.529 A is the Bohr radius and r s is the dimensionless density 
parameter (Seitz radius). We consider lowest-energy states with a given spin polarization 
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where p^ and = p—p-\ , respectively, are the (uniform) densities of spin-up and spin-down 
electrons. Including a neutralizing positive background, the total energy per electron is a 
unique function of the dimensionless parameters r s and (, 

e to t(r s ,() = t s (r s ,C) + e x (r s ,() + e c (r s ,C). (2) 
The non-interacting kinetic and exchange energies, 
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(all energies are given in units of 1 Ha = e 2 /a^ = 27.21 eV in the following), may 
be understood as the Oth- and the lst-order terms of a perturbation expansion for the 
electron-electron interaction (where r s turns out to be the expansion parameter). 

The remaining correlation energy in Eq. (J2J) appears to have the perturbation (high- 
density) expansion [TUl ITT] 



e c (r s ,C) = E MO Hr s ) + &„(C) 



n=0 



(r s < 1' 



(4) 



For the 2D electron gas (but not for the 3D one), the first coefficient vanishes, ao(C) = 0. 
Consequently, the second-order (n = 0) term is e^(() = &o(C)j representing the high- 
density (r s — > 0) limit of e c (r s ,£). It can be split into an exchange ("26") and a ring- 
diagram ("2r") term [Hj, 



4 2) (C) = ef)+ef)(C). 
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The exchange term has only equal-spins contributions, e< 26 > =e$ } (0 + ^(0, 

given 



by the 5 CT1(T2 term of Eq. (14) in Ref. (we choose the k x axis in the direction of q), 

l r°° dq r /• l l 



,(») (A) 



d 2 fc, 



\qe x + ki + k 2 | q + k lx + k 



(6) 



2.T 



A[K*(0,q] A[ Kct (C),<?] 

Here, q, k 1; and k 2 are dimensionless, a G {j, I}, and the domain of the 2D integrals is 



A[k, q] = { k G R 2 |k| < k, |k + ge x .| > k}, k«t(C) = 1 + s g n (°")C 
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[fto-(C) is the Fermi wave vector for spin-cr electrons in units of its value at £ = 0.] Scaling 
the integration variables by some constant k, q — nQ and k = kK, we have generally 



jH^ J d 2 kf(q,K) = * 2 jT^ / d 2 K f(Kq,nK). 
° 9 ^Ml ° A[1,Q] 

Applying this rule to the integrals in Eq. (J0J), we find jB] 

egJr(C) = [l + sgn(a)Cl J (26) . 



(8) 
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Consequently 0, the full second-order exchange term e£ 26 ) = e^?(C) + 641 (0 = 2J( 26 ) is 
^-independent. A Monte Carlo integration yields 



J (2b) = eg? (0) = (57.15 ± 0.05) mHa 



ImHa = 10~ 3 Ha). 



(10) 
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The ring-diagram term e^ r \() is the remaining part of expression (14) in Ref. [§], 
with the contributions 
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The equal-spins terms (di = cr 2 ) can be treated in the same way as the integral ©, 



1 + sgn(cr)C 



j(2r) 



= (76.69 ± 0.03) mHa. 
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The only non-trivial ^-dependence is in the opposite-spins term e^j) (£) = e cit (C)i 



(2r) 



(C) = e§?(0) [l-/(C) 
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By definition, /(0) = 0, and, since A[K|(l),g] = 0, /(l) = 1. Moreover, e§f (0) = -j( 2r ). 
When the results of a Monte Carlo evaluation of f(() at different values of ( are compared 
with the functions f a (C) = [(1 + 0° + (1 — C) a — 2]/(2 a — 2), particularly good agreement 
(specially for £ — > and £ — > 1) is found in the limit a — > 1 (Fig. la), 



/(C) = /i(C)+*/(C), 
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[Note that f a {C) a lso represents the ^-dependence of t s (a = 2) and e x (a = |) in Eq. (jSJ). 
The small deviation 5/(C) is accurately fitted by a polynomial (Fig. lb) 



5/(0 « 0.0636 C 2 - 0.1024 C 4 + 0.0389 C 6 • 



(15) 



The small minimum of 5/(0 indicated by the numerical data (dots in Fig. lb) at £ ~ 0.98 
is probably real, since a similar peculiarity is observed for the 3D electron gas (see the 
inset in Fig. 1 of Ref. [12!). 



In summary, the second-order correlation energy ef^ (() = e (2b ^ + (£) is 



,(26) 



,(2r), 



4 2j (C) = e?T(C)+2e^(0 + <|(C) = [153.38 /(C) - 192.46] mHa, 
where f(() is given by Eqs. ((H)) and (fTKj) . The spin resolution is fixed by 
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1 + x 19.54 mHa. 



(17) 



e^\Q = e c (0, is the high-density limit of the general correlation energy e c (r s ,(). 
To illustrate the relevance of this limit for finite densities (r s > 0), the present result 
can be used in the interaction-strength interpolation (ISI) of Ref. ^3]. This approach 
does not require the higher-order (n > 1) terms of the expansion (j3J) (which is expected 
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to have only a finite radius of convergence). Instead, information beyond the second 
order is taken from the low-density (strong-interaction or Wigner- crystal) limit of the 
exchange-correlation energy e xc = e x + e c (per electron), 



(r s ,C) ^7^ + ^ (rwoo). 



The coefficients [H] « —1.1061 and m | are independent of (, since any spatial 
overlap between two electrons is strongly suppressed in this limit, no matter whether their 
spins are parallel or not jT3]. The resulting ISI expression for the exchange-correlation 
energy at finite densities reads [13] 



:i + y)"._ 1 _ zln ((i + y)V ' 
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Using b^ = \ and writing e^r^) = c x (()/r s , we have explicitly [13] 
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Eq. (|T9"|) provides a simple explicit LSD, 

E L J D [ PhPl ] = J <Prp(r)e£'(r.(r),C(r)), (21) 

for treating arbitrary 2D electron systems (also finite ones such as quantum dots) by the 
Kohn-Sham Equations of DFT. In Eq. (|2*T|) . r s (r) = a^ 1 [7rp(r)] -1 / 2 and £(r) = |/>f(r) — 

P|(r)]M r )- 

In Fig. 2a, the ISI prediction e^ 5/ (r s ,^) = e 1 ^ 1 — e x for the correlation energy of the 
unpolarized uniform electron gas (( = 0) is compared with the accurate parametrization 
of the fixed-node diffusion Monte Carlo results in Ref. [Sj. e 131 differs slightly from 
the latter by up to 4%. This mild deviation might be cured by including in the ISI a 
simple model for the next-order coefficient of expansion ^ [16j . In the high-density limit 
(r s — > 0), however, where the present result is exact, the parametrization in Ref. [S] has 
for 0.7 < C < 0-95 a small positive deviation j3], shown in Fig. 2b. 
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Figure captions: 



Fig. 1. (a) Numerical results (dots) for the function /(£) of Eq. (fT3|) obtained by Monte 
Carlo integrations of expression (JTTJ) (with <J\02 =ti) & t selected values of (. The ana- 
lytical function /i(C) of Eq. (fTlj) is plotted as a dashed curve. The solid curve represents 
the accurate fit f t {() + 6f((), using Eq. (JTSJ) for <5/(C). (b) The fit (HSJ) (solid curve) 
compared to the true deviation (dots) of the Monte-Carlo-integration results from fx (£) . 

Fig. 2. The correlation energy of Ref. [3] (dotted curves) versus the present ISI results 
(solid curves). 
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Figure 1: (a) Numerical results (dots) for the function /(£) of 
Eq. (fTHj) obtained by Monte Carlo integrations of expression (JTTJ) (with 
o"iO"2 =TI) a t selected values of (. The analytical function /i(C) of 
Eq. (JT3J) is plotted as a dashed curve. The solid curve represents 
the accurate fit /i(C) + 6f((), using Eq. dHJ) for 8f{(). (b) The fit 
(fTo]) (solid curve) compared to the deviation (dots) of the Monte- 
Carlo-integration results from fi(()- 




Figure 2: The correlation energy of Ref. [5J (dotted curves) versus 
the present ISI results (solid curves). 



